function xopt = qpdantz_kmv_accelerated(H,f,A,b,xmin,maxiter,a,inv_H,TAB,rhsc,mnu,nc,ibi,ili)
%QPDANTZ solves a quadratic program using an active set method.

rhsa=a-f;

basisi=[inv_H*rhsa;
       rhsc-A*inv_H*rhsa];

% 
% mnu=length(f);
% nc=length(b);
%     
% ibi=-[1:mnu+nc]';
% ili=-ibi;
    
[basis,ib,il,iter]=qpsolver(TAB,basisi,ibi,ili,maxiter);


xopt=zeros(mnu,1);



for j=1:mnu
%    if il(j) <= 0
%       xopt(j)=xmin(j);
%    else
       xopt(j)=basis(il(j))+xmin(j);
%    end
end
